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Abstract. We present new results of point-contact Andreev-reflection (PCAR) 
spectroscopy in single-phase Mgi_ x Al x B2 single crystals with x up to 0.32. Fitting 
the conductance curves of our point contacts with the two-band Blonder-Tinkham- 
Klapwijk model allowed us to extract the gap amplitudes A CT and A ff . The gap values 
agree rather well with other PCAR results in Al-doped crystals and polycrystals up 
to x = 0.2 reported in literature, and extend them to higher Al contents. In the low- 
doping regime, however, we observed an increase in the small gap on increasing 
x (or decreasing the local critical temperature of the junctions, T^) which is not as 
clearly found in other samples. On further decreasing below 30 K, both the gaps 
decrease and, up to the highest doping level x = 0.32 and down to Tjr — 12 K, no gap 
merging is observed. A detailed analysis of the data within the two-band Eliashberg 
theory shows that this gap trend can be explained as being mainly due to the band 
filling and to an increase in the interband scattering which is necessary to account for 
the increase in A, at low Al contents [x < 0.1). We suggest to interpret the following 
decrease of A ff for < 30 K as being governed by the onset of inhomogeneity and 
disorder in the Al distribution that partly mask the intrinsic effects of doping and is 
not taken into account in standard theoretical approaches. 
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1. Introduction 

As well known, the superconductivity in MgB2 is characterized by two distinct energy 
gaps due to the presence of various bands crossing the Fermi level (generally grouped 
in two systems: the 3D 7r bands and the 2D a bands) and to the exceedingly small 
quasiparticle scattering between these bands. The values of the gaps measured in pure 
MgB 2 by many different techniques agree very well with those calculated within the 
two-band models in BCS pQ or Eliashberg approach [2j. The multi-band nature of 
MgB 2 allows explaining most of its features - in particular, the relatively high critical 
temperature and its unexpected robustness against sample quality - but also hugely 
increases the complexity of the effects that can arise when the system is in some way 
"perturbed" . 

A particularly interesting and debated point in the physics of MgB 2 is the possibility 
to attain the so-called "gap merging", i.e. the complete isotropization of the compound 
with consequent collapse of the two gaps in a single gap with BCS character pQ. 
Within the two-band Eliashberg theory, the gap merging can be attained, for example, 
by keeping all the parameters as in pure MgB 2 and only increasing the interband 
scattering rate, 7 CTVr . The two gaps approach each other asymptotically as a function 
of 7 CT7r , while the critical temperature is reduced; for sufficiently high values of this 
parameter (> 50 meV), one expects the two gaps to be virtually indistinguishable. The 
critical temperature of "isotropic" MgB 2 varies between 19 and 26 K, depending on the 
calculations [21 [1] . Actually, such values of •y a7T are not physical since it is practically 
impossible to increase the interband scattering without affecting other parameters of 
the material, namely the partial DOS of the a or tt bands. This is true for chemical 
substitutions in MgB 2 , e.g. carbon in the B site and aluminum in the Mg site [3], 
but also, unexpectedly, for neutron irradiation [I]. Generally speaking, the coexistence 
of various effects makes it difficult to experimentally single out their contributions. 
Fortunately, pair breaking from interband scattering gives rise to peculiar effects so 
that it can be separated rather easily from other sources of T c reduction: 7 CT7r suppresses 
T c and the large gap A CT , simultaneously increasing the smaller gap A n . 

In a theoretical paper by Erwin and Mazin [5], Al substitution in the Mg site was 
proposed as an effective way to increase the scattering between bands. First-principle 
calculations gave, for 2% of Al, a value of j an = 1.1 meV, which is already expected to 
have measurable effects on the critical temperature and on the gaps [5] . This made the 
Mg!_ x Al x B 2 system the most likely candidate for the attainment of the gap merging. 

In the past years, many experimental efforts have been made to test these 
predictions. Aluminum substitution was indeed one of the first successfully obtained in 
MgB 2 [3]. Al atoms in MgB 2 are almost completely ionized, exactly like the Mg atoms 
they substitute. Aluminum is thus a donor and the three electrons in its outer shell 
enter the system of bands giving rise to electron doping. The effects of the substitution 
on the lattice are rather complicated by a strong tendency to the formation of different 
phases. Early reports [3] showed the presence of two phases with A1B 2 structure and 
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different c axis in polycrystalline samples of nominal Al content between 0.1 and 0.25. 
A similar result was found in single crystals grown at ETH (Zurich) by means of a high- 
pressure, cubic-anvil technique. In this case, the precipitation of a non-superconducting 
MgAlB 4 phase was observed for x > 0.1 by High- Resolution TEM and other structural 
analysis techniques [6j. An independent, indirect confirmation to this picture came from 
point-contact Andreev-reflection (PCAR) measurements we performed in those single 
crystals, which showed an anomalous trend of the gaps A CT and A n as a function of the 
Al content [7J[_], with a crossover between two regimes around x = 0.1. For x < 0.1, the 
large gap decreased linearly with x while the small gap showed a pronounced tendency 
to increase, as theoretically expected |5j. For x > 0.1, A CT was found to saturate at 
about 4 meV, while A^ was fast suppressed becoming smaller than 1 meV at x = 0.2. 
An analysis of the data within the two-band Eliashberg theory showed that the trend 
observed for x > 0.1 could be explained by a decrease in the 7r-band superconducting 
coupling, that we argued could be related to the precipitation of the spurious phase [6]. 

These anomalies were not confirmed by successive measurements we performed 
in polycrystalline Mg 1 _ ;c Al ;c B 2 samples grown in Genova [8] that did not suffer from 
extended phase segregation even at high Al contents (x = 0.2). In fact, their growth 
technique involved a very long high-temperature reaction (100 h at 1000 °C) and no 
evidence of spurious phases was found by XRD - even though successive microprobe 
analyses (wavelength-dispersive X-ray spectroscopy, WDX) showed a small amount (4% 
of the volume at most) of a secondary phase in the form of micrometer-size islands 
embedded in a (Mg,Al)B 2 matrix [9 J that were concluded to have negligible effect on the 
superconducting properties. PCAR measurements we carried out in these polycrystals 
showed an almost linear decrease of both A CT and A^ as a function of the Al content which 
turned out to be in qualitative agreement with the findings of specific-heat measurements 
in the same samples [8] as well as with the results of PCAR measurements in crystals 
and polycrystals carried out by other groups [101 E] • I n a U these cases, no gap merging 
was observed up to x = 0.2, but its occurrence at a higher Al content (corresponding to 
T c around 12 K)pl)] was apparently suggested by the overall gap trends. 

Recently, new Mgi_a.ALj.B2 crystals have been grown at ETH that do not show phase 
segregation up to x = 0.32. In this paper we present the results of PCAR measurements 
on this new generation of single crystals, and compare them to the results of PCAR 
on polycrystals grown in Genova [8J as well as to other data in literature [TTJ], [11], [_2] . 
We will show that: i) our data extend previous results in single crystals [_0] up to the 
region of extremely high doping; ii) our data differ from most of the results in literature 
in the low-doping region, where we observe a much more marked increase in the small 
gap A,r on increasing x. Once reported as a function of the critical temperature of the 
junctions, T C A , the values of A n reach a maximum around T C A = 30 K and then start to 
decrease; iii) this trend nicely agrees with that observed by specific-heat measurements 
in high-quality polycrystals free of compositional gradients [12J; iv) our data show no 
gap merging up to x = 0.32 and down to T C A = 12 K. We will also show that the trend 
of the gaps in our single crystals can be well explained within the two-band Eliasberg 
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Figure 1. (a) The reconstructed hOl reciprocal space section of the AN394/5 sample 
(x = 0.24). (b) The reconstructed hkO reciprocal space section of the AN412/5 sample 
(x = 0.32). 



theory as being due to the band filling (which is the dominant effect of Al doping) and to 
a substantial increase in interband scattering in the low-doping region (x < 0.1). This 
intrinsic effect of Al doping explains the observed initial increase in on increasing x, 
while for x > 0.1 other phenomena, e.g. inhomogeneities in the dopant distribution - 
witnessed by a sudden increase in the superconducting transition width - may concur 
in making the gaps decrease again. In this range of doping, all the theoretical models 
that do not take into account inhomogeneity and disorder should be used with some 
caution as their predictions might not reflect actual properties of the compound. 

2. Experimental 

2.1. The samples 

Mgi_ x Al x B2 single crystals were grown by using the high-pressure cubic-anvil technique 
described elsewhere [131 [H], starting from pure B powder and a Mg-Al alloy and tuning 
time, pressure and temperature to eliminate phase segregation. As usual, the Al content 
was determined by energy dispersive X-ray analysis (EDX) and ranged between 0.02 
and 0.32. A thorough characterization of the structural, morphological, chemical and 
superconducting properties was performed, by using the techniques described in [6]. In 
particular, the quality of the crystals was checked by XRD using a X-ray diffractometer 
equipped with CCD area detector (Xcalibur PX, Oxford Diffraction), which allowed 
us to examine the whole reciprocal space (Ewald sphere) for the presence of other 
phases or crystallites with different orientation. No additional phases (impurities, twins 
or intergrowing crystals) were detected by examination of the reconstructed reciprocal 
space sections. This is clearly seen by comparing figured] of this paper to figure 5 of [6]. 

The crystals revealed MgB2 structure [15] . The structure refinement results of two 
samples with high Al content are presented in Table [TJ Because Al and Mg have almost 
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Table 1. Structure refinement and crystal data for Al-doped MgB2 samples AN394/5 
and AN412/5. 





AN394/5 




AN412/5 


Empirical formula 


/"A/Ft A 1 \ 13 






Temperature, K 
Wavelength , A / radiation 




29 j(2) 






n 7in7i / ii/fn'k'*-!/ 




Diff Tactometer 




Oxford diffraction 4- circle 

d iff r ac tome tcr ( CCD detector) 




Crystal system 








Space group 




P6/mmm 




Unit cell size (A) 


a = 3.0787(4) 




a = 3.0673(5) 


Unit cell volume (A 3 ) 


c = 3.5198(4) 




c = 3.4258(5) 


28.892(6) 




27.913(8) 










Absorption correction type 




analytical 




Crystal size (mm) 


U.00 X U.lO X U.Uo 




n n v n 00 ^ n no 

U.ol X U.ZZ X U.Uz 


Thcta range (cleg) 


o. 1 y to 00 . io 




7 RO +r\ n 

( . oy to 00. ud 


Limiting indices 


~3 < < 4 
— 5 < I < 5 




_ 4 < A < 2 
~ 2 < I <5 


Reflections collected/unique 


204/46, FUnt=0. 0364 




156/37, llmt=0. 0179 


Refinement method 




Full-matrix least-squares on .F 2 




Data/restraints/ parameters 


46/0/6 




37/0/6 


Goodness of fit on _F 2 


1.005 




1.071 


Final R indices [I > 2ct(Z)] 


Rl=0.0243, wR 2 =0.0611 




Rl=0.0162, wR 2 =0.0409 


R indices (all data) 


R!=0.0324, wR 2 =0.0620 




Rl =0.0174, wR 2 =0.0411 


A Prnax , Ap min (e/A 3 ) 


0.354 and -0.350 




0.222 and -0.247 




Fractional atomic coordinates 


and atomic displacement parameters (A 2 ) 




B 




x=l/3; y=2/3; z=l/2 




(Mg,Al) 




x=0; y^0; z=0 




Un, B 


0.010(1) 




0.007(1) 


U33, B 


0.014(1) 




0.010(1) 


U 12 , B 


0.005(1) 




0.003(1) 


Un, (Mg,Al) 


0.009(1) 




0.008(1) 


U33, (Mg,Al) 


0.012(1) 




0.008(1) 


U 12 , (Mg,Al) 


0.005(1) 




0.004(1) 



the same amount of electrons (12 and 13 respectively) the refinement was performed 
without Al and the position of Mg was considered to be occupied with both atoms. 
The presented refinement results of two highly-doped samples and their reconstructed 
reciprocal space sections show that the crystals chosen for subsequent PCAR analysis 
are of a high quality and satisfy the requirements for single crystals. On the basis of 
these data, we can exclude the influence on measured properties of some possible factors 
like polycrystallinity, additional phases, strong disorder etc. 

Figure [2] shows the relation between the Al content and the bulk T c given by DC 
susceptibility in our single crystals (black circles). The critical temperature was defined 
here as the abscissa of the intersection between the y = axis and the linear fit of the 
susceptibility vs. T in the region of the transition. The transition width 5T C (defined 
as T 10 % — T 90 %) increases on increasing the Al content, from 0.77 K in the crystal with 
x = 0.02 up to about 9 K in the most doped sample (x = 0.32). The large values of ST C 
in heavily doped samples can be related to the local inhomogeneity in the dopant content 
[16] , and to the smallness of the coherence length £ that allows inhomogeneities on 
a scale of the order of £ to be resolved p2]. A contribution from the simple disorder 
consequent to Al doping may be present as well, but other sources of broadening are 
either negligible or excluded by the single-crystal nature of our samples. For instance, 
due to the small DC field used in the magnetization measurements (2 to 5 Oe), effects 
related to the magnetic field are certainly small. It will be clear in the following that 
the inhomogeneous distribution of Al in the crystals does not invalidate the results of 
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Figure 2. Bulk critical temperature T c measured as a function of the aluminum 
content x. • : Data obtained in our single crystals by DC susceptibility [B]. A: Data 
from [10], obtained in single crystals by specific-heat measurements. □: Data from DC 
susceptibility in long-annealed polycrystals [18] . 



local measurements such as PCAR spectroscopy, provided that the gap amplitudes are 
reported as a function of the local critical temperature in the region of the contact and 
not as a function of the (average) Al content or of the bulk T c . 

2.2. Point-contact Andreev-reflection measurements 

Before starting with the point-contact measurements, we cleaned the crystals and etched 
their surface by dipping them into a solution of 1% HC1 in dry ethanol. After 2-5 
minutes, we rinsed the crystals in pure ethanol and dried them with nitrogen. The point- 
contact measurements were performed by using the "soft" technique described elsewhere 
[T9] . Instead of pressing a metallic tip against the sample as in standard PCAR, we made 
the contact by using as a counterelectrode a small spot of Ag conductive paint. This 
pressure-less technique can be used also on brittle samples and, in thin single crystals, 
allows injecting the current (mainly) along the ah planes [19], so as to measure both the 
a and 7r-band gap at the same time [20] . The diameter of the Ag-paint spot is typically 
< 50 /iin that, however, does not correspond to the actual size of the point contacts. 
As a matter of fact, parallel microjunctions are very likely to form between the crystal 
surface and the Ag particles in the paint within the macroscopic contact area, so that the 
measured I— V characteristics and conductance curves should be regarded as an average 
over a certain region in direct space. Usually, the potential barrier at the N/S interface 
is rather low so that the contacts are in the Andreev-reflection regime. Otherwise, the 
characteristics of the contact and its normal-state resistance can be tuned by using short 
voltage or current pulses [4, 19J. The formation or modification of contacts with the help 
of electric pulses is well known in standard electrotechnics [21]. It was also used, already 
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Figure 3. Temperature dependence of two conductance curves measured in point 
contacts on crystals with x = 0.24 (a) and x — 0.32 (b). A subset of the complete 
series is shown for clarity. The temperatures are indicated in the labels. The bottom 
curve in each panel is the normal-state conductance and the temperature at which it 
is reached is defined as the Andreev critical temperature Tf. 

in the Seventies, to create point contacts for phonon spectroscopy in normal metals [22] 
or high-quality Josephson contacts between two superconductors [23] . During a voltage 
pulse (in our case of the order of 1 V for some milliseconds), the contact region can 
be heated well above the bath temperature. This phenomenon was shown to give rise 
to local annealing in heavily neutron- irradiated MgB2 [I], but in the present case this 
drawback can be ruled out because the crystals are already well annealed [6] and no 
enhancement of the local T c above the bulk value has ever been observed. 

Figure [3] shows two examples of raw conductance curves measured in the crystals 
with x = 0.24 and x = 0.32 as a function of temperature. A subset of the complete 
series is shown for clarity. The values of the normal-state resistance of the junctions 
are Rn = 120 Q and = 178 Q, respectively. Such high values of Rjy are necessary 
to fulfill the conditions for ballistic transport through the junction [24], because of the 
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shortening of the mean free path due to Al doping. Let's suppose for a moment that a 
single contact is established between the crystal and a Ag grain in the paint; then the 
whole resistance is due to this single contact and, using the residual resistivity of the 
most-doped crystal (p — 5 /ifi-cm), the relationship between p and the mean free path 
£ reported in [25], and the Sharvin formula [21], one can evaluate the contact radius a 
from the contact resistance. If — 180 Q the result is a ~ 11 A that is indeed much 
smaller than the mean free path £ ~ 100 A. This ensures ballistic transport even if a 
single contact is established between sample and counterelectrode. If several contacts 
are present, as it probably happens in our case, each of them has higher resistance than 
the parallel as a whole and is thus certainly ballistic. 

The temperature at which the Andreev-reflection features disappear and the 
normal-state conductance is recovered will be in the following referred to as the local 
critical temperature of the contact, or the 'Andreev critical temperature" T C A . The 
values of T C A are reported in figure [3] for both the contacts. In doped samples, because 
of the local inhomogeneity in the Al content, different contacts on the same sample can 
provide different gap amplitudes and different T^. All the values of T C A are included 
between the onset and the completion of the magnetic transition, so that in samples 
with a wide superconducting transition (especially the most doped ones) values of T C A 
substantially smaller than the bulk T c can be obtained. For these reasons, is more 
appropriate than the bulk T c to describe the properties of the contact. 

3. Results 

Figure H] shows the conductance curves G(V)=dI /dV measured in a6-plane contacts 
in single crystals with different Al contents x from up to 0.32. All the curves are 
normalized, i.e. divided by the normal-state conductance (measured at T = T c or in a 
magnetic field H = H c2 ). The rather small amplitude of the normalized conductance 
curves is related to the "soft" point-contact technique we use. In particular, it is due to 
some very thin (smaller than the coherence length £) impurity layer on the surface of 
the Ag grains in the paint, which gives rise to inelastic scattering at the interface. As 
shown in (26] , this effect does not affect the measured gap values and can be accounted 
for by simply inserting an extrinsic broadening in the BTK fit of the conductance curves 
(see Appendix A). 

To obtain the gap values, we fitted the experimental curves with the two-band 
BTK model in which the conductance G of the junction is a weighed sum of G a and G n : 
G = w n G n + (1 — w^Go- [20J. G a and G n depend on the relevant gap amplitude (A CT 
or A n ), on the effective potential barrier parameter (Z a or Z n ), and on a broadening 
parameter (r CT or r„-), in conformity to the conventional BTK model [27] modified 
by including the effect of the quasiparticles lifetime [281 121] (for further details, see 
Appendix A). The two-band BTK best-fitting curves are indicated in figure H] as solid 
lines. 

In the low-doping regime (up to x = 0.15) there is no point to compare the two- 
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Figure 4. Normalized low-temperature conductance curves in crystals with different 
Al contents from x = to x = 0.32. The Andreev critical temperature and the 
normal-state resistance Rn oi each contact are also indicated. Circles: experimental 
data. Solid (dashed) lines: best-fitting curves obtained within the two-band (single 
band) BTK model. 

band BTK fit with the single-band one; the latter is always of poor quality and does not 
reproduce the position of the peaks and the width of the Andreev features. However, in 
the most doped samples (where, in principle, a gap merging could occur) this comparison 
cannot be omitted. As a matter of fact, dashed lines in the two bottom panels of the same 
figure represent the single-band BTK curves that best fit the experimental conductance 
curves. As discussed in Appendix A, the two-band fit is always preferable (even in the 
x = 0.32 case, where the two theoretical curves are very similar to each other) on the 
basis of a statistical Fisher F-test. This would lead to the conclusion that two gaps are 
always present in the Al-doped crystals, even at the highest doping content. 

However, the most reliable test for the actual presence of two gaps consists in 
studying the magnetic-field dependence of the conductance curves, owing to the faster 
suppression of the 7r-band gap by the magnetic field [30j. In pure MgB 2 , this technique 
allowed us to separate the partial a and 7r-band conductances and to fit each of them 
with a standard, three-parameter BTK model [19], [30]. In doped samples, the complete 
separation is not always possible but, if two gaps are present, an outward shift of 
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Figure 5. Magnetic- field dependence of the conductance curves measured on the 
crystals with x = 0.24 (a) and x = 0.32 (b). The magnetic field was applied parallel 
to the ab plane. Vertical lines indicate the shift in the position of the peaks. 



the conductance maxima is observed at a certain magnetic field, when the a-band 
conductance becomes dominant. Figure [5] reports the magnetic- field dependence of 
the conductance curves of the two contacts in the most Al-doped crystals (x = 0.24 
and x = 0.32) whose temperature dependencies have been already shown in figure [3j 
Vertical lines indicate the maximum shift of the conductance peaks, that witnesses the 
presence of two gaps (rather close to each other) and justifies the two-band BTK fit 
reported in the two bottom panels of figure HI 

The gap amplitudes A CT and A n given by the two-band fit of the conductance curves 
of our point contacts in Al-doped single crystals (of which figure 0] showed a subset) are 
reported as a function of the Al content x in figure E](a). The vertical spread of data 
for each doping content gives an idea of the variation in the local gap values in different 
contacts on the same crystal. The trend of the gaps A a and is however clear: The 
large gap monotonically decreases on increasing x in the whole doping range, while the 
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Figure 6. (a) Energy gap amplitudes A a (O ) and A„. (• ) measured by PCAR in 
single crystals, as a function of the Al content determined by EDX. Error bars indicate 
the uncertainty on the gap values for a given curve, (b) The same values as a function 
of the critical temperature of the contact T^. The gaps (obtained by PCAR) from refs. 
[5] (squares), [TU] (triangles) and [TT] (down triangles) are also reported for comparison. 
The dashed lines are only guides to the eye, while the straight grey line in (b) indicates 
the BCS A vs. T c dependence. 



small gap first slightly increases - reaching a maximum of 3.1 meV at x — 0.08 - and 
then starts to decrease. For x > 0.15, the slope of the two curves is apparently the 
same. 

The vertical dispersion of data can be partly removed by plotting the gaps as a 
function of the critical temperature of the contacts, T C A , as in figure 0(b). The same 
figures also report other PCAR data from literature, obtained by us in polycrystals [8] 
(squares) and by other groups in single crystals [10] (up triangles) and polycrystals [TT] 
(down triangles). Recent STM measurements of the 7r-band gap in crystals grown at 
ETH substantially agree with these data [31 J . 

In the high-doping region (T C A < 25 K) all the data sets agree very well with one 
another independent of the nature of the samples. In particular, our data nicely extend 
the curves previously obtained in single crystals by Klein et al. down to very low critical 
temperatures (well below the theoretical limit for gap merging in the hypothesis of pure 
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interband scattering [TJ, 32J). In the low-doping region, a substantially common trend is 
observed for the large gap A a , although our data (both in crystals and polycrystals) are 
a little higher than those reported in literature [HlJ E] • One may wonder whether this 
shift is related to the "soft" PCAR technique. This is not the case, since in pure MgB 2 
we obtained gap values in excellent agreement with those shown here (at the highest 
T C A ) also by standard PCAR measurements with Au or Pt tips [33] . Whatever its origin, 
this shift makes our values for A CT and A^ fall above and below the BCS value (straight 
line in figure [6b), respectively, while in many data reported in literature A a is very close 
to the BCS line or even falls below it. A similar effect has been recently observed in 
heavily neutron-irradiated MgB2, but in that case it can be somehow related to the high 
level of disorder, on the basis of analogous evidences for conventional superconductors. 
Here, this anomaly is very difficult to justify theoretically - and is certainly beyond 
current models for two-gap superconductivity. 

As far as the small gap A„. is concerned, our data show that, in the low-doping 
region, it increases on increasing the Al content, reaching a maximum around = 30 
K. This effect is quite definitely assessed by the measurements we expressly carried out 
in crystals with x = 0.02 and 0.034. Such a tendency is also present in the specific-heat 
data from [12] and, although much smaller, in the data from [10] and [11] as well as 
in the results of specific-heat measurements in polycrystals from Genova [8j [M]. The 
tendency of the small gap to increase was recognized in [18J as an intrinsic effect of Al 
doping, more evident in samples produced via a long reaction at high temperature so 
as to reduce the strain and the inhomogeneity in the Al content. The much greater 
increase in A n at low doping contents in the single crystals grown at ETH (similar to 
that of samples of the "B" series in [H]), might then be related to the absence of lattice 
strain due to compositional gradients, also witnessed by the rather sharp transition 
of the crystals in this range of doping levels. If this picture is valid, the increasing 
inhomogeneity of the crystals on increasing x might be responsible for the fact that, for 
T C A < 25 K, the small gap of these single crystals returns on the same curve described 
by the values of A„- in the other samples. 

4. Discussion 

We tried to interpret the trend of the g function of the local critical temperature 

within the two-band Eliashberg theory. A model for the effect of Al doping on the gaps 
of MgB2 was given by Kortus et al. [32] who solved the two-band Eliashberg equations 
scaling the Eliashberg functions by the change of the DOS alone. The same approach 
was used in [35] to explain the x dependence of the critical temperature in Al-doped 
MgB2- This model has no free parameters as long as one takes the interband scattering 
rate to remain negligible as it is in pure MgB 2 (see Appendix B for the explicit Eliashberg 
equations and the details of the model). However, even in this case it proved sufficient 
to qualitatively explain the previous experimental data for the gaps in Al-doped MgB 2 
[341 E] as a function of the critical temperature, thus indicating that the changes in the 
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Figure 7. Main panel: Energy gap amplitudes in our Al-doped MgE>2 single 
crystals (• ,0 ) as a function of the Andreev critical temperature Tjr . Gap 
amplitudes from other PCAR measurements are reported for comparison 
(squares [8], up triangles [TO] and down triangles [H]). Lines indicate the 
predictions from Eliashberg theory, when only the DOSs are changed according 

to the band filling due to electron doping ( ), and when a proper increase 

in the interband scattering at intermediate doping levels is also included in 
the model (— — ). The straight grey line represents the values of a BCS gap. 
Inset: doping dependence of the gaps in our single crystals, compared to the 
Eliashberg curves (DOS + interband scattering). 



a- and 7r-band DOS are by far the dominant effect of Al doping. 

Dashed lines in figure [7] represent indeed the T c A -dependencies of the gaps A CT and 
A^ we calculated within the two-band Eliashberg theory following the aforementioned 
approach [32J ES] and using the DOS N CT (E^) and N^E^) from first-principle 
calculations [36]. It is clearly seen that the dashed lines already reproduce rather well 
the data we obtained by PCAR in polycrystals [8] (squares) and those reported in Refs. 
[TO] and [UJ. As far as the PCAR data in our single crystals are concerned (circles in 
figure UJ), it is clear that the large gap would be perfectly compatible with the DOS 
scaling alone (dashed lines), while the small gap is definitely not. As a matter of fact, 
the initial marked increase in on decreasing T A requires some additional ingredient 
in the model. According to the discussion of section [TJ and to the predictions of [5], 
this immediately suggests an increase in the interband scattering parameter 7^ due to 
the Al doping (incidentally, let's recall that the intraband scattering parameters, which 
certainly increase with x, are however ineffective in changing the gap(s), according to 
Anderson's theorem. Strictly speaking, this is only true as long as the doping content 
is small so that the perturbative description of the doped compound is possible). By 
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taking 7 CT7r as the only adjustable parametejj] and using the calculated DOS (as we did 
for the dashed lines in figure [7]), we were able to reproduce both the critical temperature 
and the gap values of our single crystals. The resulting curves for A CT and A^ are shown 
as solid lines in figure [7J 

To obtain these curves, a non-monotonic dependence of 7 CT7r on T C A was necessary. 
7 CT7r initially increases (almost quadratically) on decreasing T C A , reaching a maximum 
7™ = 2.6 meV at T C A ~ 29 K (i.e., x ~ 0.1), then decreases linearly to finally saturate 
to 0.22 meV at T C A = 15 K. At T C A = 33.5 K (that means x ~ 0.02), the value of the 
interband scattering parameter used for the fit (7^ = 1.0 meV) is perfectly compatible 
with the theoretical predictions (7^ = 1.1 meV)[5]. 

If the initial increase in 7 CT7r on increasing the Al content is easily explained in terms 
of out-of-plane distortions of the B sublattice [5] , the problem arises of explaining the 
decrease in j an at higher Al contents that is necessary to reproduce the observed gap 
trend in single crystals, which shows no gap merging down to T C A values as small as 
12 K. Giving a definitive answer to this problem certainly requires further theoretical 
and experimental investigations. Here, a simple interpretation can be anticipated. As 
discussed above, the comparison of our experimental data to those of [12] and [T8] 
suggests that the initial increase in A n , and the corresponding increase in interband 
scattering, are intrinsic effects of Al doping in MgB 2 , as theoretically predicted. At the 
end of section [31 we have also proposed to interpret the decrease in A n for T C A < 30 K 
(x > 0.1) as being due to the onset of inhomogeneity in the Al content. As a matter 
of fact, just in the same region the slope of 5T C as a function of T c (being ST C the 
width of the superconducting magnetic transition) suddenly increases. The simplified 
Eliashberg model we have been using (and which is described in detail in [35] as well as 
in Appendix B) is unsuited to take into account these effects. As a matter of fact, it is a 
mean-field model that treats the doped MgB 2 as a perturbation of the pure compound 
so that: i) the Anderson's theorem holds; ii) the mathematical expression for the matrix 
elements of the Coulomb pseudopotential and of the coupling constant derived for pure 
MgB2 holds as well. This may not be true any longer when the doping concentration is 
too large: in this case the model itself probably fails and a different, non-perturbative 
description should be used. In this sense, the decrease in 7^ for T C A < 30 K necessary 
to fit the experimental A^ values may not reflect an actual property of the system. 
In other words, this decrease might be necessary to mimic the effects of lattice stress 
and inhomogeneity (possibly at a nanometric scale) in the local Al content that, at the 
present moment, are not explicitly included in the model. Incidentally, this situation is 
somehow similar to what we observed in heavily neutron- irradiated MgB2 [4]. 

t Actually, a very small change in the prefactor of the Coulomb pseudopotential, /io, was necessary as 
well to consistently reproduce the correct values. However, /io varies from 0.031 at T C A = 39 K to 
a minimum of 0.027 reached when = 30 K. 
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5. Conclusions 

In conclusion, we performed a large number of point-contact Andreev-reflection 
measurements on segregation-free, state-of-the-art single crystals of Mg 1 _ a ,Al ;c B2 
extending previous PCAR results up to x = 0.32 and down to local critical temperatures 
T C A ~ 12 K. The local critical temperature of each contact was directly obtained from the 
temperature dependency of the conductance, while the gap amplitudes were determined 
by a two-band BTK fit of its low-temperature bias dependency. No merging of the 
energy gaps has been observed down to the lowest T C A and the persistence of two gaps 
at the highest values of x has been confirmed by studying the conductance curves in 
the presence of suitable magnetic fields. When compared to the theoretical results 
obtained in the framework of the two-band Eliashberg theory, the experimental gaps 
show that the main effect of Al doping is to fill up the bands, thus changing the DOSs 
at the Fermi level just as expected by first-principles calculations. Nevertheless, at 
an intermediate aluminum content, corresponding to T C A values between 18 K and 35 
K, the A 7r (T c j4 ) curve of our single crystals shows clear deviations from the theoretical 
behaviour expected for pure band filling. Within the two-band Eliashberg model already 
used to describe pure MgB2, these deviations can only be reproduced by introducing 
a proper amount (up to 2.6 meV) of interband scattering 7 CT7r . However, / y a7T must 
be decreased again on further increasing x to account for the experimental behaviour 
of and the absence of gap merging - at least down to T C A = 12 K. We propose 
to interpret this trend as resulting from two competing phenomena: i) the increase 
in interband scattering, intrinsic to Al doping but clearly observable only when the 
lattice stress due to compositional gradients is eliminated, and ii) the gradual onset of 
inhomogeneity when the Al content is increased above a certain threshold (x ~ 0.1). 
At high doping contents, inhomogeneity probably dominates so that any theoretical 
description that does not take it into account becomes less and less satisfactory. To 
overcome this problem, which is likely to occur in any doped MgB 2 -based system, a 
detailed experimental knowledge of the kind of disorder would be required, as well as 
new and specific theoretical approaches. 
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Appendix A. The two-band BTK fit 

In this Appendix we will give some additional details about the two-band BTK fit and 
the fitting procedure. In the theoretical model we used, the normalized conductance 
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through the point contact is expressed by: 

G = w n G n + (1 - w v )G a (A.l) 
Each conductance is expressed in the form 



Z^/2 a sA E ^ ( t ) ) cos(0)d0 
£^/2°X*(^) cos (0) d 



Gm = j -=^p^ /,rT.yr r ( A - 2 ) 



where i = cr, 7r and: 

cos(0) 
cos(0) 2 + Z, 



= . rl (a. 3) 



|1 + (o"JV,i(0) - lJi^JT 

The functions Fi(E) are given by 

^7 



and contain the broadening parameters r CT and T v as imaginary parts of the energy 
Also note that these parameters are independent of the gap values Aj. This model is a 
two-band generalization of the formulation by Kashiwaya et al. [38J that reduces to the 
simplest BTK formulation if one takes Tj = and = instead of integrating over the 
angle as in equation IA.21 

In the modified BTK model (281 129] T is a measure of the intrinsic lifetime 
broadening. In our case, T a and T n account for both the intrinsic lifetime broadening and 
other effects - related to the experimental technique and thus extrinsic - that smooth 
the curves [19J. The most probable origin of the additional broadening is inelastic 
quasiparticle scattering in the vicinity of the contact, i.e. in a degraded layer covering the 
Ag grains of the paint. As recently shown [26] this scattering can be simply accounted 
for by increasing the broadening parameter T in the BTK model. Z a and Z n depend 
on the potential barrier at the interface. The weight of the 7r-band conductance, w n , is 
taken as an adjustable parameter as well. 

The range of variability of the seven fitting parameters is actually limited by some 
physical constraints. For example, in pure MgB 2 must vary between 0.66 (for pure 
a6-plane tunneling) and 0.99 (for c-axis tunneling) |20j. Owing to the non-perfect 
directionality of PCAR, we always took 0.68 < w n < 0.75 for a6-plane current injection 
|19j . Moreover, w n and the barrier parameters Z a and Z n must be independent of 
temperature and magnetic field so that one is forced to keep them constant in fitting 
the whole T or if-dependence of a conductance curve. In principle, IV and T n should be 
smaller than A CT and A^, respectively, and they must increase on increasing the applied 
magnetic field [30J. 

Finally, the uncertainty on A,,- and A^ for a given curve can be defined as the 
maximum range of gap values that allow a good fit of the curve, when the other 
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Figure Al. (a): normalized experimental conductance curve of a point contact in a 

crystal with x — 0.32 (O ) compared to the single-band ( ) and two-band ( ) 

BTK fit. (b) The values of A (• ), T (A) and Z (□) for the it band, that allow fitting 
the same curve, in a series of different fits, (c) The same as in (b) but for the a band. 
The variations in the gaps are qualitatively indicated by grey strips. 



parameters are changed too. To define which is the best fit, we minimized the sum 
of squared residuals (SSR). This corresponds to minimizing the chi-square but does not 
require an estimation of the uncertainty on the conductance for each point, which can 
vary from curve to curve and is often difficult to estimate. We then allowed a variation 
of the SSR of the order of 100% and determined the corresponding range of parameters. 
In the best cases, as in pure MgB 2 (top panel of figure Hj), A^ is directly related to 
the position of the peaks in the conductance curves, and A^ to the shoulders on the 
sides. In this situation, the uncertainty on the gap values is usually rather small (of 
the order of 0.3 meV). In doped samples, clear conductance peaks are still present but 
no structure at V > V pea k is directly visible that would reveal the presence of a second 
gap. Moreover, in heavily doped crystals, the two gaps can be so close to each other 
that the conductance peaks can occur at some intermediate energy. In all these cases, 
a statistical test (the Fisher F-test) can be used to determine whether the single-band 
or two-band fit is preferable. In practice, one first determines the best fitting curves 
within the two models. If the SSR of the single-band fit is smaller, this fit is certainly 
preferred. However, if the two-band fit gives a smaller SSR value, the F-test allows 
testing whether, within a fixed confidence level (usually 5 %), the improvement in the 
fit is not only due to the increase in the number of parameters from 3 to 7. In the two 
bottom panels of figure HI the best-fitting single-band BTK curve is represented by a 
dashed line. In the most ambiguous case (x = 0.32), the SSR in the range [—10, 10] mV 
(excluding the noisy regions of the curve) is 5.0 • 10 -4 for the two-band fit and 1.1 • 10~ 3 
for the single-band one. The F-test shows that, for any level of confidence, the two-band 
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fit is preferable. Another example, always in the high-doping limit (x = 0.32, T C A = 12.3 
K), is reported in figure IAT1 Here the best SSR in the range [—10, 10] mV is 1.75 • 10 -3 
in the single-band case and 9.67 • 10~ 4 in the two-band case. The F-test again shows 
that the two-band fit is better for any level of confidence. The panels (b) and (c) of 
the same figure give an idea of the spread of the fitting parameters in different fits of 
the same curve; the range of variation of the gaps is indicated by the grey strips. It 
is clear in figure IAl( b) that is comparable or equal to A^. This is a drawback of 
our experimental technique, due to the aforementioned extrinsic broadening. It must be 
said, however, that contrary to a widespread belief, the gap structures are clearly visible 
in the conductance curves (and a reliable gap measure can be extracted from their fit 
[T9] ) even if A = T, provided that Z is not too small, as in our case. This can be easily 
shown by calculating the conductance curves within the modified BTK model. 

Appendix B: Two-band Eliashberg equations 

Let us start from the generalization of the Eliashberg theory [39] for systems with 
two bands, which has been already used with success to study the MgB 2 system 
[20| El HQ]. To obtain the gaps and the critical temperature within the s-wave, two- 
band Eliashberg model one has to solve four coupled integral equations for the gaps 
Aj(iu;„) and the renormalization functions Zi(iu n ), where i = a, n is the band index 
and u n are the Matsubara frequencies. We included in the equations the non-magnetic 
impurity scattering rates in the Born approximation, 7^. 

u n Zi(ko n ) = u n + 7rT^ Aij(iuj n - \u m )N z (xu m ) + 

m,j 

+ !><^(iwn) (B-l) 

3 

Z i (i^ n )A l (i^ n ) = 7rTy^[A i:? -(ia; n - iu m ) - //*-(o; c )] • 

m,j 

■ 6{\w c \ - u m )N A (iu m ) +J2^ N L( iuJ n) (B.2) 

3 

where 9 is the Heaviside function, uo c is a cut-off energy and Aij(iu n — iuo m ) = 
J +o ° duja^{uj)F(uj)/[(uj n -uj m ) 2 +uj 2 }, N 3 A (iu m ) = Aj(h} m )/ yju^ + A?(iu m ), N J z (iu) m ) = 

Um/yju>m + A|(iu; m ). 

The solution of the Eliashberg equations requires the following input: i) four (but 
only three independent) electron-phonon spectral functions o^Au)F{oj)\ ii) four (but 
only three independent) elements of the Coulomb pseudopotential matrix fi*(u c ); iii) 
four (but only three independent) impurity scattering rates 7^. 

The four spectral functions a 2 j(u)F(u) were calculated for pure MgB 2 in ref. [ID]. 
For simplicity, we will assume here that the shape of the aj 'Auj)F(uj, x) functions does 
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not change with x, and we will only scale their amplitude with the electron-phonon 
coupling constants A^-: 



Xtj(x) : <4 i F(u,,x = 0) 



Xij(x = 0) 



where 



(B-3) 



(B.4) 



N J N (E F ,x = 0) 

As far as the Coulomb pseudopotential is concerned, we used the expression 
calculated for pure MgB 2 |JT], though including the dependence of the densities of 
states at the Fermi level N l N (E F ,x) on the doping content x: 



A* (*^)7rij (^)7r7r 
2.23 



N°(E F ,x) N°(E F ,x) 



1 



2.48 



N«(E F ,x) NUE F ,x) 



(B.5) 



As for the scattering rates, let us remind here that, according to Anderson's 
theorem, the intraband scattering parameters ja have no effect on either T c or the 
gaps so that they can be dropped. We are thus left only with the interband scattering 
rates 7 CT7r and 7^, which are however related through the equation 

7 = 7 (B6) 

This is the reason why we can choose 7 CT7r as the only adjustable parameter. 
Moreover, it can be shown that only the interband scattering can make the small 
gap increase while both T c and the large gap decrease (on increasing x), in the way 
experimentally observed. 

It is worth clarifying the relationship between the broadening parameters r a and 
in the BTK model and the scattering parameters 7 CT(J , 7^, r y a7T in the Eliashberg 
theory. The intrinsic BTK linewidth for a given band, e.g. IV, takes into account all the 
scattering channels and is thus proportional to 7^ + 7^0-. A direct relationship between 
the intrinsic T n and 7 CT7r cannot be established unless an independent determination of 
the intraband scattering rates is obtained. Furthermore, in our case the values of T n 
extracted from the fit contain an "extrinsic" term (related to inelastic scattering in the 
vicinity of the interface) in addition to the intrinsic linewidth. This prevents a direct 
connection between the Tj parameters of the BTK model and the 7^ of the Eliashberg 
theory. 



Appendix C: A test of consistency 



With reference to figure [?l it is worth noting that we are fitting with the Eliashberg 
theory - including interband scattering - data points that were obtained through a BTK 
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fit of experimental curves (and thus without taking into account the possible effect of the 
interband scattering on the curves themselves, which is theoretically expected not to be 
negligible [12] )• This approximation is hardly avoidable since a fit of the conductance 
curves within the two-band Eliashberg theory would be a very complex (if possible) 
task. For the same reason, a direct proof of the reasonability of the BTK approach - 
that would be obtained by directly comparing the results of the two fitting procedures, 
BTK and Eliashberg - cannot be obtained. However, one can try to demonstrate that 
a theoretical curve calculated within the two-band Eliashberg theory (with given values 
of the gaps and a non-zero interband scattering rate) can be fitted with the BTK model, 
and that the gap amplitudes resulting from the fit are consistent with the original ones. 

Let us refer for clarity to the point where the 7r-band gap is maximum (see figure [7]), 
which is also the most critical one. We first calculated the gap functions A a (uj) 
and A n (uj) {u being the energy) within the Eliashberg theory, with the parameters 
corresponding to the two curves in figure [7J (dashed and solid lines), that means 
in particular with 7 CT7r =0 and 7 CT7r =2.40 meV, respectively. Then, we calculated the 
corresponding Andreev-reflection conductance curves at 4.2 K with no additional 
smearing, and using the experimental values of Z a and Z v [38]. The two curves are 
reported as open circles in figure [CTl fa) and (b), respectively. 

Finally, we fitted these curves with the two-band BTK model and compared the 
values of the gaps given by the fit with those used to generate the curves. Note that 
the amplitude of the theoretical curves in figure IC1I (> 1.4) is greater than that of 
the experimental ones, which are further smeared by extrinsic broadening factors (e.g. 
related to the specific measurement technique we use [131 26J). This makes the test even 
stricter. The best-fitting BTK curves are shown in figure ICl( a) and (b) as solid lines. 
It is clearly seen that both the gaps are rather well re-obtained, with an error which is 
of the same order of magnitude as the experimental uncertainty. This shows that the 
trend of the small gap A n obtained by the BTK fit of the conductance curves is not due 
to an artifact introduced by the BTK fit itself. In other words, the enhancement of A n 
at intermediate T C A values is a real effect that would be obtained as well by fitting the 
experimental conductance curves with the more appropriate Eliashberg theory. 
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Figure CI. (a) Symbols: the theoretical (normalized) Andreev-reflection 
conductance curve generated within Eliashberg theory by using the parameters 
corresponding to the point at = 30 K on the dashed curves of figure [71 The 
values of the parameters are indicated on the left. Solid line: best fit of the 
curve with the two-band BTK model. The fitting parameters are indicated on 
the right, (b) Same as in (a), but with the parameters corresponding to the 
point at = 30 K on the solid curves of figure [7J 
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